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CLUSTER EXPANSION IN THE CANONICAL ENSEMBLE 



ELENA PULVIRENTI AND DIMITRIOS TSAGKAROGIANNIS 

Abstract. We consider a system of particles confined in a box A C M'' interacting via a 
tempered and stable pair potential. We prove the validity of the cluster expansion for the 
canonical partition function in the high temperature - low density regime. The convergence is 
uniform in the volume and in the thermodynamic limit it reproduces Mayer's virial expansion 
providing an alternative and more direct derivation which avoids the deep combinatorial 
issues present in the original proof. 



1. Introduction 

The quantitative prediction of macroscopic properties of matter through its microscopic 
structure has been a main challenge for statistical mechanics. In this direction a very impor- 
tant theoretical as well as practical contribution is the work of J. E. and M. G. Mayer [13] 
(see also Ursell [19]) in the theory of non-ideal gases where they derive a full series expansion 
correcting the equation of state for the ideal gas (p = /cTp, where p is the pressure of the 
system at temperature T with density p, k being the Boltzmann constant) to all orders in p 
obtaining the famous virial expansion. Convergence of this series has been addressed later 
(see [9] and [16]). The main idea in [13], (see also [12]) is to describe all possible interactions 
among the particles of the non-ideal gas by representing them as linear graphs, which has 
later led to a main tool in statistical mechanics, namely the cluster expansion method. The 
thermodynamic pressure is computed as the infinite volume limit of the logarithm of the 
grand canonical partition function, which is however a function of the activity of the system. 
To get an expansion with respect to the thermodynamic density, one needs to further ex- 
press the latter in a power series of the activity, invert it and replace it in the equation for 
the pressure. This gives the virial expansion after using (for the inversion) some interesting 
combinatorial properties of enumeration of graphs which finally express the coefficients of 
the virial expansion to be sums over only "irreducible graphs". 

This road which leads to an expansion of the free energy versus the density is evidently 
not the straight one! The direct and natural way is to take the density p instead of the 
activity as order parameter and correspondingly to work in the canonical rather than in 
the grand canonical ensemble, which however rests on the possibility to cluster expand the 
canonical partition function. We came to this problem from other directions, see the end 
of the introduction, and found, to our surprise, that the problem is not only solvable but 
easy. It fits in fact beautifully in the theory of cluster expansion for abstract polymer 
models, as developed in all details by many authors after the pioneering work of [7], [6]. The 
exponentiation procedure in this theory produces a lot of diagrams which are not present in 
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the Mayer expansion and which therefore must vanish in the thermodynamic hmit. As we 
shall see in Section 5.1 the origin of such cancellations is closely related to the basic property 
of the cluster expansion (from which the expansion takes its name) namely that the only 
chains of graphs (clusters) which survive in the expansion are made of "incompatible" graphs. 
The validity of the cluster expansion for the canonical ensemble opens the way to attack 
several other problems (which was actually our initial motivation) such as the finite volume 
corrections to the free energy, the radius of convergence of the expansion in powers of the 
density (rather than the activity) for both the general model and the particularly interesting 
case of hard spheres, the construction of coarse-grained Hamiltonians via multi-canonical 
constraints as required by the Lebowitz-Penrose Theorem [10] for Kac interactions with 
applications to the LMP model [8] and its variants. We hope to further address these issues 
in subsequent papers. 



2. The model and the result 

We consider a configuration q = {gi, . . . , qj^j} of particles in a box A := (— C M'^ 

(for some L > 0) which interact with a stable and tempered pair potential : M*^ ^ M, i.e., 
there exists B > such that: 

5^ V{q,-q,)>-BN, (2.1) 

l<i<j<Af 

for all A^ and all qi, ...,qN and the integral 

C(/3) := [ le-^^^-'^-llrfg (2.2) 

is convergent for some /3 > (and hence for all /3 > 0). Since in this paper we are interested 
in the infinite volume limit of the free energy density, we can assume periodic boundary 
conditions since it is a general result (see e.g. [16] and [5]) that the thermodynamic limit is 
independent of the choice of the boundary conditions. This particular choice in the present 
paper is not essential, periodic boundary conditions are used in order to obtain translation 
invariance in some cases (see e.g. Lemmas 4.1 and 5.2). Furthermore, our result remains 
valid with other boundary conditions by slightly changing the proof and we hope to address 
it in a subsequent work where we will consider the finite volume corrections as well. 
We obtain the periodic boundary conditions by covering M.'^ with boxes A and adding all 
interactions. Let 

V^'''{q^, q,) := ^(^^ " 'i^ + '^^)' (2-3) 

where apart from assuming stability and temperedness, we further need to guarantee conver- 
gence of the above sum by imposing a condition on the decay properties of V . A potential 
V is called lower regular if there exists a decreasing function ip : ]R_|_ — IR+ such that 
V{x) > —iIj{\x\) for all x G M"' and Jg°° ?/'(s)s^~^(is < oo. Then V will be called regular if it 



CANONICAL ENSEMBLE 3 

is lower regular and there exists some ry < oc such that V{x) < ipHxl) whenever > ry 
and this is our extra assumption (see also the discussion in [5]). 

The canonical partition function of the system with periodic boundary conditions is given 
by 

^/ dg,...dg^e-^K-i^\ (2.4) 
iV ! Jj^N 

where is the energy of the system with periodic boundary conditions given by 

^r(q):= E ^"^'■(?-?.)- (2-5) 

l<i<j<N 

Given p > we define the thermodynamic free energy by 

f^ip) := hm /;3,A(iV), where f^^^iN) := --^ logZ^X^- (2-6) 

The main result in this paper is: 

Theorem 2.1. If pC{l3) is small enough (see (4.6)^ then 

1 I A 1^ 

log Z^;i^ = log ^ + m E (2.7) 



where the coefficients FN^A{n), n > 1, satisfy 

|i^iv,AHI <Ce-^", (2.8) 
for constants C, c > that are independent of N and A. Moreover, 

FNAin) = ^— Piv.|A|(n)5^,A(n), (2.9) 

where -P/v,|A|(n) and Bp^\[n) are defined in (5.2) and which in the thermodynamic limit give 

lim |y^|(n) = p" and \\m Br kin) = (2.10) 

7V,|A|-s>oo,7V/|Al=p ' A-s>oo 

for all n > 1, where fin are the irreducible coefficients of Mayer: 

f^-'-=^ J2 [ n ie-^'^^'^^-'^'^ -^)dq2---dqn+i, Qi = (2.11) 

and Bn+1 is the set of2-connected (or irreducible) graphs g on (n+1) vertices (i.e., connected 
graphs which by removing one vertex and all related edges remain connected) and E{g) is the 
set of edges of the graph g. 

Note the unfortunate coincidence of notation between the inverse temperature /3 and the 
irreducible coefficients which however we keep in agreement with the literature. To 
prove Theorem 2.1 we first establish in Section 3 the set-up of the cluster expansion for 
the canonical partition function in the context of the abstract polymer model following 
[1], [14] and [7]. Note that we could also work with more general formulations such as in 
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[15], [3] or [17]. In Section 4 we prove the convergence condition and as a corollary of the 
cluster expansion theorem we prove (2.8). The discussion of the thermodynamic limit is 
left for Section 5 where using the splitting in (2.9) we prove (2.10) by investigating the 
cancellations that lead to the "irreducible" coefficients given in (2.11). This latter fact 
is the result of some special structure (that we call "product structure") which may occur 
already in the abstract formulation of the polymer model and we address it in full generality 
in Subsection 5.1. Last, we also give an appendix with the main ideas of the original virial 
expansion and connect it to our approach. There is a number of interesting combinatorial 
issues that arise, many of which have been extensively studied in graph theory. With no 
intention of being exhaustive we just refer to [11] which is inspired by the virial expansion. 



3. Cluster expansion for the canonical partition function 

We view the canonical partition function Z^^^ ^ as a perturbation around the ideal case, 
hence normalizing the measure by multiplying and dividing by |A| in (2.4) we write 

7per yideal yint 

where 



yper yideal yint ("i ^\ 



For Z'^l^\j^ we use the idea of Mayer in [13] which consists of developing e~^^^ '^'^^ in the 
following way 

e-/^-r(.)= n (1 + A.)=E n (3-3) 

l<i<j<7V gGSiv {i,j}GE{g) 

where by we denote the set of all simple graphs (undirected with no self-loops and no 
multiple edges) on vertices (the labels of the particles), where a graph is the pair g = 
(y{g), E{g)) where V{g) C {1, . . . , A^} is the set of vertices and E{g) := {{z, j}, z, j G V{g)} 
is the set of (non-oriented) edges, and 

e-^^''^^fe-*)-l. (3.4) 

Given a graph g E Qn say that it is connected if any two vertices of g are connected 
by following its edges. We define the support of a graph (denoted by supp g) to be the set 
of its vertices. Two graphs g,g' are compatible (denoted by (7 ~ g') if suppf? fl supprjf' = 
(i.e., they are not connected); otherwise we call them incompatible {'y^). We also denote 
by \g\ the cardinality of supp (7. With these definitions we see that any graph g E Qn '^^ 
equivalent to the pairwise compatible (non-ordered) collection of its connected components, 
i.e. g = {gi, . . . , (7fc}~ for some k. Hence, 

k k 
{gi,...,gk}^ *=1 {Vi,...,Vfe}^ i=l 
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where 

UiV) := Yl Ca(^?), Ug) ■■= I n m n (3.6) 

Here Cy denotes the set of all connected graphs with support in V C {1, . . . , A^}. We will also 
be using the notation C„ for the set of all connected graphs on n vertices. Both expressions 
in (3.5) are in the form of the abstract polymer model which we specify next. 

An abstract polymer model (F, Gr, w) consists of (i) a set of polymers F := {71, ...,7|r|}, (ii) 
a binary symmetric relation ~ of compatibility between the polymers (i.e., on F x F) which 
is recorded into the compatibility graph G-p (the graph with vertex set F and with an edge 
between two polymers 7i,7j if and only if they are an incompatible pair) and (iii) a weight 
function w : F — )■ C. Then, we have the following formal relation which will become rigorous 
by Theorem 3.1 below (see [7], [1] and [14]): 



where 



or equivalently ([1],[2]) 



Cl 



Wi-ir^^\ (3.8) 



GcGi 



Cl 



n <9^('ri)a;(7i)---a^(7")a;(7„) 



a;(7)=0 



(3.9) 



The sum in (3.7) is over the set X of all multi-indeces / : F — )■ {0, 1, . . .} with supp / := {7 G 
F : 1(7) > 0}, uj^ = (7) is the graph with '^■yi^suppi ^i'^) vertices induced 
from ^supp/ C Gr by replacing each vertex 7 by the complete graph on 1(7) vertices. 
Furthermore, the sum in (3.8) is over all connected subgraphs G of Qi spanning the whole 
set of vertices of Qi and /! = n7esupp/ -^(7)'- Note that if / is supported on a compatible 
collection of polymers 7 then c/ = 0. 

We state the general theorem as in [1], [14] but in a slightly different form. Let 

L = L(S) = =up (zMi^l = zMl^, (3.10) 

where 6 will be small and therefore L = 1 + 0(5). The optimal bound for the convergence 
radius is beyond the scope of the present paper, however, we hope to come back to this issue, 
also in the particular case of hard spheres, in a subsequent work. 

Theorem 3.1. [Cluster Expansion] Assume that there are two non-negative functions a, c : 
F — )■ M such that for any 7 G F, |a;(7)|e'^*-'*''' < 5 holds, for some 6 > small. Moreover, 
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assume that for any polymer 7' 

^|a;(7)|e'^W+'=W<ia(7'), (3.11) 

where L is given in (3.10). Then, for any polymer 7' G F the following bound holds 

J2 |cX|e^^e-pp^^(^)^(^) < L|a;(7')|e"(^')+"(^'), (3.12) 

/:/(7')>l 

where Cj are given in (3.9). 

Proof. Apply the proof in [1] for the activities uj{'j)e^^"'\ □ 

In view of (3.5) we can represent the partition function Z^^j^^ both as a polymer model on 

graphs with (C(A^), Gc, Ca) where we explicit the dependence on N by denoting by C{N) the 
set of all labeled connected graphs on up to vertices and as a polymer model on vertices, 
with (V(iV), Gv, Ca) where V{N) := {V : V C {1, . . . , A^}}. 



4. Convergence of the cluster expansion, proof of (2.8) 

In this section we check the convergence condition of Theorem 3.1. We work in the case of 
vertices which corresponds to the abstract polymer model (V(A^), Gy, Ca) where V(A^) : = 
{V : y C {1, . . . , A^}}. Then, as a corollary of Theorem 3.1 we prove (2.8). 

Lemma 4.1. There exist two positive functions a,c : V(A^) — t- R such that for any V G 
V{N): 

|CA(^)|e"(^)<5 (4.1) 
holds for some S > small. Moreover, for any V G V(A^) 

E ICA(F)|e'^(^)+'=(^)<ia(0, (4.2) 



L 

V: VrpV 



where L is given in (3.10). 



Proof. We choose aiV) = a\V\ and c(y) = c\V\, for some a, c > to be chosen later and let 
a := a + c. To bound |Ca(^)| we use a version of the tree-graph inequality (proved in this 
form in [15] and [16]) which states that for a stable and tempered potential, we have the 
following bound: 

|E n /M|<e^^""E n (4-3) 

g£Cu {i,j}eE[g) rer„ {i,j}eE(T) 
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where 7^ C is the set of trees with n vertices. Then 



|a(V)|e«l''1 < e'™-")!"!/ nil E H I/mI 

! jgy I I TGTJvl {i,j}(iE(T) 

< e(^^--)"5: / W---W n (4-4) 
Given a rooted tree T let us call (ii, ji), {12,32), {in-i,jn-i) its edges. We have: 

n-l 



/ w ■ ■ ■ w n = / ■ ■ ■ n i/^-iJ 

JA" |A| |A| |A| JAn 

^ / dq,, [ dy2-- - [ dy^fl\e-^y-^y^^-l\ 
'^1 J A J A J A r_o 



~ A 



< 



iAr 



fc=2 
n-l 



|A| ^ 



lAI" 



C^a(/3)' 



(note that the choice of V^'^^ makes Ca(/3) independent of x) where we considered qi^ as the 
root and we used the change of variables: 

Vk = qi, - qj„, yk = 2, ...,n. (4.5) 

We choose pC{f3) such that: 

^^2/3B+a+l^(^^ < 1. (4.6) 

Then, since the number of all trees in Tn is from (4.4) we obtain (recalling that 

N/\A\=py. 

n-2 1 

lCA(V^)|e'^l^l < ^e(^^^+'^)"CA(/3)"-^ < ^pCA(/3)e^(^'^^+'^), (4.7) 

by using the bound 2 < n < N and the fact that pe^^^+''CA(/3) < 1. The latter is 
true considering that inequality (4.6) still holds with Ca(/3) for A large enough, since 
limA^ooCA(/3) = C{/3). Then defining 6 := |pC(/3)e2(2/3s+«), (4.1) is satisfied. Consid- 
ering a = 1 we have (for any fixed i): 

recalling that L = 1 + 0(5) and where S' := pe^'^^+''+^C{(3). Since 7^ V'} C Uievi'^ ^ i} 
we get (4.2) and conclude the proof of the lemma. □ 

The way we chose to present the cluster expansion as well as its convergence can by no means 
give the best radius of convergence. Our goal was merely to obtain (giving up the seek for 
the best radius) the consequence of the cluster expansion theorem, given in (3.12), which we 
use in order to establish (2.7). Nevertheless, our condition is comparable with the ones in 
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the literature (see [9], equation (3.15), and in [16], Theorem 4.3.2) and we will clarify these 
issues in a subsequent work. Moreover, for the particular case of the hard spheres we can 
obtain the improved radius as in [4] but for the density p rather than the activity, since we 
are working directly with the canonical partition function. 

Proof of (2.8). After proving the convergence condition it follows directly from Theorem 3.1 
that for all V E V{N) and by choosing c{V) := c\V\ and a{V) := \V\ 

E |c/Cl|e'="^ll<L|CA(0|e"l^'l ||/||:= ^(^)l^l' (^.9) 

r.i{V')>i vesupp/ 

where we remind that a = 1 + c. Thus, exponentiating the partition function (3.5) using 
(3.7), we obtain (2.7): 

i^iog^S,^ = ^E^^cl = ^ E _ E , ^'^c1 = ^E^^.aH' (4-10) 



.. ,N} r.supp lB{i} ' ' ' ' n>l 

where 



7:supp79{l} 

|/|=n+l 



U y 

yGsupp/ 



(4.11) 



Note that |/| < ||/||. The function -F/v.A(n) is uniformly bounded for all A^, A as well as 
absolutely summable over n, namely from (4.9) with V"' = {1} we get: 



'N,A 



'n) < 



Y E |c/Cl|e^" < e-^"Le-, (4.12) 



n . _ 

/:supp/3{l} 

|/|=n+l 



which concludes the proof of (2.8). □ 

Note that the sum in (4.11) is infinite. As it will be apparent in the next section it is more 
convenient to work with a truncated version of F^^/^ln). Thus, for M > ri + 1 we define 

F^'4n):=^ E ^^Cl, (4.13) 

i':supp/3{l} 
|/|=n+l, ||/||<M 

which is uniformly bounded for all A^, A and absolutely summable over n using the same 
argument as in (4.12). Note also that now the sum is finite. Furthermore, for all n > 1 

-f(n+Af) 

sup |F^a(^) - FNA{n) \ < sup y |c/Cl|e5("+*^) < e-i("+*'^)Le", (4.14) 

N,A ' ' Af,A n+1 r:.s^x 

/:supp/9{l} 

|/|=n+l, ||/||>M 

since inside the sum we have: n < M < \\I\\. Hence, we can next work with the truncated 
function Fff^Jji) and then pass simultaneously to both limits: the thermodynamic and 
M oo. 
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5. The thermodynamic limit, proof of (2.10) 

The sum in the definition of -F^^^(?7,) (4.13) does not depend on the labels of the extra n 
particles (we have already chosen label 1). Hence, 

^U^^ = V^li^ n^) E c,Ci = ^^PN,Un)Bl[Anl (5.1) 

\ / r- II ni^ H4 



I: \\I\\<M 
supp/={l,...,r!,+l} 



where 

PN,Un) := ^^"\ -i^""^ and 5,-f^(n) := ^ (5-2) 

supp /={l,...,n+l} 

While obviously P/v,|A|('^) p", for Bp^j^{n) we proceed as follows: we have restricted 
ourselves to a system of only ra + 1 particles, therefore we can view the sum in B^^pjji) as 
the cluster sum in {y{n + 1),Gv,Ca) with the extra condition ||/|| < M. While for the 
convergence it was more convenient to work with the space (V(iV), Gy, Ca); the cancellations 
that we discuss in this section occur already at the more elementary structure of labeled 
connected graphs (i.e., on (C(n + 1),Gc,Ca) whose multi-indeces we denote by / : C{n + 
1) — !■ {0, 1, . . .}) so since ||7||<jv/ '^/Ca = Ylil\\i\\<M ^i'^h (coming from the same partition 
function) we rewrite B^fJ^ri) as follows: 

= ^ E <=/c.; = E E '/cl. (5.3) 

/:||/||<M £/eC„+i /:||/||<A/ 

Note that looking at the support of some cluster / we see a connected graph (with possible 
superposition of some graphs g ii I{g) >2, oi with some common parts of a graph if they are 
incompatible in more than one points). Thus, in (5.3) we first put in evidence the resulting 
graph g G Cn+i and then sum over all / that correspond to it. The main idea of what follows 
is to investigate the order in |A| of the products Ca- We start with a definition: 

Definition 5.1. Given a connected undirected graph a vertex is said to he an articulation 
point if removing it and all edges incident to it the graph results in a non- connected graph. 

One of the main ingredients of our proof is the fact that if {g' : g' G supp /} are connected at 
only articulation points of the fixed g G C„+i then the activities Ca factorize. To emphasize 
this fact we state it as a lemma: 

Lemma 5.2. Given g G and an incompatible collection G such that Ug'^cd' = 9; if for 
all incompatible pairs g',g" we have that supp fl supp (7" = {i} is a singleton and moreover, 
the vertex i is an articulation point of g then 

Ug) = n (5-4) 

g'SG 
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for all finite A. 



Proof. By the definition of an articulation point tliere is a. g' & G wliicli lias only one 
articulation point. We integrate the coordinates in g' keeping fixed the one that corresponds 
to the articulation point. In principle the result should depend on the position of the 
articulation point, but this is ruled out by the special choice of boundary conditions (as we 
stressed in Section 2) which makes such integrals translation invariant. Hence, the result of 
the integration gives Ca(5'')- Then, again by the definition of an articulation point, this can 
continue until integrating all g' E G. □ 

Coming back to (5.3) the power of |A| is n — X]g'6supp/(lfi''l ~ ^)^i.9')- Moreover, since it is 
always true that n + 1 < ^^/gsupp/dfi''! — 1) + 1 (by the fact that all g' G suppJ should be 
incompatible, i.e., have at least one common label) it is implied that non-negligible terms 
(in the limit A — >■ oo) should satisfy 

I{g') = 1, ^g' G suppJ, and (5.5) 

n + l= + (5-6) 

g'Gsupp I 

(In the representation of hyper-graphs by bi-partite graphs, these two requirements corre- 
spond to graphs without loops.) Again, to emphasize we state it as a lemma 

Lemma 5.3. Given I, if there is A G suppJ such that (\{Ugi^Ag') 7^ Ylg'eA^^id') then the 
cluster that corresponds to I is vanishing at the thermodynamic limit. 

Proof. Suppose that there is such a A G I, i.e.. 



Us' 

g'eA 



g'eA 



1. 



(5.7) 



Then, (5.7) implies that 



n + 1 = 


U 


< 




+ 






g'esupp/ 




g'eA 




g'Gsupp I\A 



g'eA 



g'Ssupp I\A 



(5.8) 

hence (5.6) cannot be true and it has to vanish by taking the limit. □ 

Proof of (2.10). Motivated by these facts we proceed with (5.3) by splitting the sum over / 
into the following two complementary cases: 

(1) (factorization property holds) (\{Ug^Ag) = YigeACA{9) fo^' all A C suppJ, 

(2) there exists A C suppJ such that (\{Ug(zAg) 7^ Y[geACA{g)- 

Lemma 5.3 states that in the limit |A| — )■ 00 the terms of Case (2) vanish (and they are 
finitely many due to the truncation). For the case (1), if g G Cn+i \ ^n+i then all terms in 
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B^\{n) are exactly canceling each other. This is a property of the combinatorial coefficients 
multiplying products of activities of polymers which have the factorization property, just as 
in the cluster expansion where compatible collections (carrying the factorization property) 
do not contribute in the sum. In the next subsection we investigate this property that we 
call "product structure" (see Definition 5.1) and since it can be a general property of the 
abstract polymer model we prove it in its general form in Lemma 5.5 and then we explicit 
it in our case in Corollary 5.6. 

Thus, the only remaining terms are those of Case 1 for G Bn+i which due to condition 
(5.5) give non-vanishing contribution only for I{g) = 1 (and those with higher multiplicity 
are only finitely many, due to the truncation). Hence, for all n > 1 we can pass to the 
thermodynamic limit (note that (4.14) permits to exchange sum and limits) and obtain: 

I Al" 

hm B^,A{n) = hm lim 5*'^(n) = lim ^ V Uig) = Pn. (5.9) 

where /3„ is defined in (2.11), concluding (2.10) and the proof of Theorem 2.1 (pending the 
proof of Corollary 5.6). □ 

5.1. The "product structure". Coming back to the general formulation (r,Gr,cj), for 
any F' C F we define the set of all incompatible sequences that can be constructed out of 
elements of F' by: 

F'^ ■={A: Ac F', incompatible}. (5.10) 

Recall that every single element 7 G F' is considered incompatible, hence the singleton {7} 
is an element in F^). 

Definition 5.4. Given (F,(Gr,w), we say that F^ C F has a "product structure" if 

• there exists a one-to-one function : F^ — > F, with 0({7}) = 7, if 1 ^ T^; 

• for any A G F^, we have the factorization 

a;(0(A)) = n ^(7')- (5.11) 

7'GA 

We also define the range of by 

R(0):={0(A),VAGFt}. (5.12) 

We are interested in all multi-indeces / such that every 7 G supp/ is the image of some 
A G F^, i.e., supp / C R(0) and this is the content of the following lemma. 

Lemma 5.5. IfT^cT has a product structure then in the expansion (3.7) we have: 

cjLo^ = cico^. (5.13) 

/:supp7cR{0) /:supp/={7'},7'er'' 
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Proof. Given r** C F with product structure, using (5.12), we define 

k 

Z*{T^) := ^r,{a;(7)=o,V7^RW} = X] ri'^*^^^) 

{7lv,7fc}~ «=1 
7i6R(0) 

k 

= E n^('^(^^))= n(i+'^(7'))- (5.14) 

The first equahty of (5.14) is due to the fact that (p is one-to-one, i.e., for any 7^ e R(0) 
there is a unique G with 0(Ai) = 7^. Then using the factorization property, i.e., 
uj{(j){Ai)) = Yl'y'eAi'^i'y')^ second equahty is due to the fact that n7'Gr''(-'- +'^(7')) = 
X^Acrf- Yl-y'eA'^i^')' "where the latter sum is over subsets A (compatible or incompatible) of 
the set of vertices in F^, hence it can be uniquely decomposed into k compatible components 
A = {Ai, . . . , Ak}^ with Ai e Fj, for all i. 

Then if we exponentiate the last expression of (5.14) we obtain the right hand side of (5.13), 
while if we exponentiate -^r,{w(7)EE0,V7^R((/>)} (by first exponentiating using (3.9) and then 
evaluating) we obtain the left hand side of (5.13). □ 

Now we explicit the result of Lemma 5.5 for the case of (5.3), which corresponds to the 
particular polymer model (C(n+1), Gc, Ca)- Any g G Cn+i\Bn+i can be uniquely decomposed 
into 2-connected components, say g = uf^ibi where bi G B{n + 1), Vz = l,...,k (recall 
the difference of notation between i3„,+i and B{n + 1) with the latter being the set of 2- 
connected graphs on up to n+1 vertices). Then, according to Definition 5.4, any such element 
{bi, . . . , bk} has product structure with 4>{A) := UbeAb for all A C {61, ... , bk} incompatible 
and it also satisfies (see Lemma 5.2) the factorization property by construction. Thus, 

Corollary 5.6. Given g G Cn+i \ Bn+i, then 

E = 0, (5.15) 

/: |j/||<A/ 
U rQ'=Q 

g tsupp 1^ ^ 

for all truncations M and where * is to remind that the sum is over I who are supported on 
polymers g which have the factorization property, as in Case 1 in the proof of (2.10). 

Proof. Suppose that {61, ... , bk} is the unique decomposition into 2-connected components 
of g and hence has product structure. Then applying Lemma 5.5 for F^ := {61, . . . , bk} we 
obtain that all coefficients of ChipiY^ ■ ■ ■ ChipkY* (with -|- . . . -|- n^l^fcl < M) are zero 

except those with k = 1, which however correspond to the case g = bi E Bn+i and are not 
considered since we assumed g G C^+i \ Bn+i- □ 
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6. Appendix: Mayer's virial expansion 

The approach introduced by Mayer, see [13], is to work in the grand canonical ensemble 
which is a measure on both the number of particles N and the configuration q 

G^,,,A(iV; dq) := ^Ar^z"" e-^^'^^^^dq, . . . dq^, (6.1) 
where a(2^) is the grand canonical partition function given by 

H/3,a(z) := Yl ^'^^/3,A,JV (6.2) 

Af>0 

and z is the activity of the system. The thermodynamic pressure is defined as the infinite 
volume limit of the logarithm of the grand canonical partition function: 

P/3(^) := , lim P(3,a(,z), where /3p/3,a(z) = log S^,a(2;). (6.3) 
The idea in [13] consists of developing e"^'^**-''-' in the following way 

,-^H.i.)= n (i+/m)=e n (6-4) 

l<i<i<7V gggjv {i,j}eE{g) 

where by Qn we denote all simple graphs on vertices, E{g) is the set of edges of a graph 
g E Qn and 

A, :=e-^^('^-*)-l. (6.5) 
Then the grand canonical partition function becomes 

S/3,a(2;) = ^ ifA(^), vfhere WA{g) := JJ /^jH^^*' ^^■^'> 

where by \g\ we denote the cardinality of the graph g and we define it to be the number of 
vertices. Using the fact that the weight w\{g) is multiplicative on disconnected components 
a general theorem on enumeration of graphs gives (see e.g. [18] where it is stated as "The 
first Mayer Theorem") 

^a(^) = exp <^ ^ ^ 5^ WA{g) \ , (6.7) 

where C„ is the set of connected graphs on n vertices. This is the predecessor of the Cluster 
Expansion Theorem 3.1! Then defining 



\A\n. 



14 ELENA PULVIRENTI AND DIMITRIOS TSAGKAROGIANNIS 

(which is normahzed in the volume and hence it has a limit 6„ := lim|A|^oo ^r!,(A)), equation 
(6.3) gives 

^' ' n>l ^ n>l 

In the thermodynamic hmit the canonical free energy is the Legendre-Fenchel transform of 
the pressure, namely 

/3//3(p) = sup{plog2; - (3pfs{z)} = p\ogz{p) - (3pp{z{p)), (6.10) 

z 

where z{p) is given by the inversion of the relation p = zp'p{z). Note that this is also 
equivalent to first defining the finite volume density by 

p^{z):=¥.G,,AN]=zp'p^^{z) (6.11) 

and then pass to the limit |A| — oo. In [18] this inversion is referred as "The second Mayer 
Theorem" and it is again a result on enumerating connected and 2-connected graphs, where 
the latter means all graphs which cannot be reduced to connected graphs by removing a 
point and all related edges. Under again the assumption that w^ig) is multiplicative (see in 
our case Lemma 5.2), we have that 



dz \ ^ n\ ^ ] \ dp \ Tn\ 

\n>l g£C„ / k \m>2 geBr, 



where Bm is the set of 2-connected graphs on m vertices. Note that this is the combinatorial 
counterpart of our discussion in Section 5. From (6.12) we have that 

p = zp'p{z) ^ z{p) = pe~^->2^'"-i'''"~\ (6.13) 

where 

hm -1- J2 ^a(^)- (6-14) 
|AKoo A m! ^-^ 

Plugging into pf^{z{p)) we obtain the famous virial expansion: 

f3p,ip)=p-J2^^^P"'^'- (6-15) 

m>l 

Overall, (6.10), gives 



I {Pi^ogp- 1)-Y: ;;^Prnp-A. (6.16) 
m>l ) 
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